home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / POLDIV.DEM < prev    next >
Text File  |  1991-05-01  |  1KB  |  38 lines
  1. PROGRAM d5r3(input,output);
  2. (* driver for routine POLDIV *)
  3. (* (x-1)**5/(x+1)**3 *)
  4. CONST
  5.    n=6;
  6.    nv=4;
  7. TYPE
  8.    glnarray = ARRAY [1..n] OF real;
  9.    glnvarray = ARRAY [1..nv] OF real;
  10. VAR
  11.    i : integer;
  12.    u,q,r : glnarray;
  13.    v : glnvarray;
  14.  
  15. (*$I MODFILE.PAS *)
  16. (*$I POLDIV.PAS *)
  17.  
  18. BEGIN
  19.    u[1] := -1.0; u[2] := 5.0; u[3] := -10.0;
  20.    u[4] := 10.0; u[5] := -5.0; u[6] := 1.0;
  21.    v[1] := 1.0; v[2] := 3.0; v[3] := 3.0; v[4] := 1.0;
  22.    poldiv(u,n,v,nv,q,r);
  23.    writeln;
  24.    writeln('x^0':10,'x^1':10,'x^2':10,'x^3':10,'x^4':10,'x^5':10);
  25.    writeln;
  26.    writeln('quotient polynomial coefficients:');
  27.    FOR i := 1 to 6 DO write(q[i]:10:2);
  28.    writeln;
  29.    writeln('expected quotient coefficients:');
  30.    writeln(31.0:10:2,-8.0:10:2,1.0:10:2,0.0:10:2,0.0:10:2,0.0:10:2);
  31.    writeln;
  32.    writeln('remainder polynomial coefficients:');
  33.    FOR i := 1 to 4 DO write(r[i]:10:2);
  34.    writeln;
  35.    writeln('expected remainder coefficients:');
  36.    writeln(-32.0:10:2,-80.0:10:2,-80.0:10:2,0.0:10:2)
  37. END.
  38.